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Abstract—This article attempts to organize and review 
the variety of different types of hierarchical vector bases 
appropriate for use in solving the common equations of 
electromagnetic fields on (possibly curved) triangular or 
quadrilateral domains. Tables are used to illustrate 16 
different basis families for triangular domains and 7 
families for quadrilateral domains. The choice of a 
suitable basis may be dictated by the relative linear 
independence of the set and the ease with which functions 
can be used on meshes containing both cell shapes. We 
also consider suitability for use with surface integral 
equation applications. 


Index Terms—Basis Functions, Boundary Elements, 
Electromagnetic Fields, Method of Moments (MoM), 
Vector Finite Elements. 


I. INTRODUCTION 


For more than 35 years, vector basis functions have found 
application in electromagnetics for representing electric and 
magnetic fields within numerical solutions of differential 
equations and representing surface current densities within 
solutions of the electric and magnetic field integral equations 
(EFIE and MFIE, respectively). The present review will 
consider both applications but maintain a primary focus on 
two-dimensional basis functions that can be used for surface 
integral equations. Early basis functions were low degree 
polynomials that were often interpolatory, with coefficients 
that represent specific current components at interpolation 
points within the mesh. Since interpolatory functions have 
the disadvantage that they are not well-suited for adaptive p- 
refinement procedures, a number of hierarchical basis sets 
have also been introduced. In the following, we review 
many of the hierarchical vector basis functions proposed 
during the past two decades to show their similarities and 
identify possible advantages of different families. 

Most of the proposed basis functions are of the curl- 
conforming variety, but we note that in two dimensions these 
functions are easily converted into divergence-conforming 


functions by a 90 degree rotation. Previous research 
suggests that the optimal basis for the EFIE is a divergence- 
conforming set that spans the mixed-order Nedelec spaces 
[1-3], while the optimal basis type for the MFIE may not be 
divergence-conforming and apparently should span 
polynomial-complete spaces [3-5]. Thus, in the following 
we will consider both types: families that span complete 
polynomial vector spaces and those that span the mixed- 
order Nedelec spaces (“reduced gradient” spaces for curl- 
conforming or “reduced curl” spaces for divergence- 
conforming) originally defined in [1]. In the case of 
triangular cells, it is not always immediately obvious which 
domain space is properly spanned by a particular family; 
many papers refer to “Nedelec spaces” without clearly 
distinguishing between the mixed-order spaces of [1] and the 
polynomial-complete spaces of [6]. (The term “mixed” is 
also misleading in that it is sometimes used to denote that 
different polynomial degrees are used in different directions 
or sometimes to denote representations that work when more 
than one cell shape — say, triangular and quadrilateral — 
arise in a mesh.) 


The families under consideration are those from references 
[7-24] for triangular domains and [11, 17, 20, 23-28] for 
quadrilateral domains. For future reference, we will tabulate 
specific functions up to order 3 from many of the proposed 
families in a common notation. The procedure for mapping 
these bases to curved cells (triangular or quadrilateral, 
divergence-conforming or curl-conforming) has been 
described in the literature [29] and will not be discussed 
here. 


П. A PERSONAL NOTE 


The lead author first investigated hierarchical vector bases 
for triangles in 1997, in connection with a sponsored 
research project involving adaptive p-refinement techniques. 
At that time, the papers of Jon Webb and his students were 
the principal sources of information on hierarchical bases 
[30-32]. For example, the 1993 paper by Webb and 
Forghani [30] presented basis functions up to polynomial 
degree 2. While this set did not contain a subspace that 
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exactly spanned the reduced gradient Nedelec space of 
degree 1.5, something that the author desired to do, it was 
relatively easy to build upon Webb and Forghani’s idea to 
create a similar hierarchical set that did properly span the 
reduced gradient spaces. That set was published some years 
later as an Appendix to [15]. 


Webb went on to publish a comprehensive set of curl- 
conforming hierarchical bases for triangles and tetrahedra in 
1999 [10]. It appears that that set was primarily intended to 
span polynomial complete spaces, but it had the nice 
attribute that the functions were organized into gradient and 
rotational subspaces (an approach employed below for many 
of the basis families), which is helpful when dealing with the 
reduced-gradient spaces. The Webb functions did not satisfy 
the conditions for spanning the reduced-gradient Nedelec 
spaces of order 2.5 or higher. At about the same time, 
Andersen and Volakis published a hierarchical set that 
spanned the reduced-gradient spaces of Nedelec, but did not 
span polynomial complete spaces [8-9]. In the years since, a 
number of other hierarchical basis families have been 
proposed, some of which span one type of space, and others 
that span both types. 


The principal limitation of many of these hierarchical 
families is the degradation of their linear independence as the 
order increases. In 2008, the present authors decided to 
collaborate on developing hierarchical basis sets that 
exhibited better linear independence than some of the 
preceding sets. That led to the publication of [21, 28, 33-34] 
and the recent text [29]. In the course of our investigations, 
we organized other available families of bases into a 
common notation and implemented many of them to test 
their performance. Since the papers written on this subject 
tend to employ different notations, we felt that a review 
article that presented the different families in a common 
notation, and pointed out their similarities and differences, 
might be useful — hence the present communication. 


III. MIXED-ORDER NEDELEC SPACES 


In our work we employ a strict definition of Nelelec’s 
mixed-order spaces of functions, namely the definitions 
presented in [1]. On triangles, these spaces discard all 
degrees of freedom associated with the gradient (or curl) of a 
higher degree polynomial (for curl conforming or divergence 
conforming, respectively). On quadrilateral cells, some 
additional degrees of freedom are added to make the 
resulting bases rotationally symmetric. We employ a half- 
integer index to denote the reduced-gradient spaces. Table 1 
shows the number of independent functions residing over 
triangular and quadrilateral cells for a particular order. 


IV. TRIANGLE BASIS FAMILIES 


Table 2 summarizes the literature on hierarchical vector 
basis families for triangular cells. These may be curl- 
conforming or divergence conforming in the original papers, 
which are easily converted into the other by a 90 degree 
rotation. Functions described on tetrahedrons are easily 
projected onto triangles. In our summary tables, we have 


converted all families into the curl conforming form for 
easier comparison. 

























































































Table 1 
Number of independent basis functions 
overlapping a 2D cell as a function of order, for 
polynomial-complete vector spaces and reduced 
Nedelec spaces 
triangle 
order complete reduced 
0 2 
0.5 3 
1 6 
1.5 8 
2 12 
2.5 15 
3 20 
3.5 24 
4 30 
4.5 35 
quadrilateral 
order complete reduced 
0 2 
0.5 4 
1 8 
1.5 12 
2 16 
2.5 24 
3 28 
3.5 40 
4 44 
4.5 60 

















Functions that span polynomial complete spaces 


J. P. Webb and his colleagues were the first to propose 
hierarchical vector basis families for triangular cells [7, 10, 
30-32]. Reference [7] presented a general expression for 
divergence-conforming hierarchical basis functions, which 
we have converted into curl-conforming functions in Table 
3. In the tables to follow, simplex coordinates (u, v, w) are 
employed, where 0 < u < 1, 0 < v < 1, and w = 1 — u – у. 
Note also that the functions labeled “edge” functions in 
Table 3 and other tables that follow straddle two cells and 
must be adjusted to maintain tangential-vector continuity 
across that boundary (for the curl conforming case), while 
functions labeled “face” are local to a cell. 
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Table 2 


Some of the available families of vector bases on 
triangular cells, with indication of the spaces they 
are designed to span. 





reference complete reduced 
[7], 1997 X 
[8-9], 1998-9 X 
[10], 1999 
[11], 2001 
[12], 2001 
[13], 2001 
[14], 2003 
[15], 2004 
[16-17], 2005-6 
[18], 2006 
[19], 2007 
[20], 2007 
[21], 2011 
[22], 2011 
[23], 2013 
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The functions of Table 3 span polynomial complete 
spaces. In reference [10], Webb presented a different family 
of curl-conforming bases that also span polynomial complete 
spaces, with the functions organized into subgroups of 
rotational and gradient functions (Table 4). Webb also 
considered the orthogonalization of the basis functions to 
improve the matrix conditioning associated with their use, 
and the functions that take a different form (the higher order 
rotational functions) are presented for degree 3 in Table 5. 
Comparing the basis functions of [7] and [10], we observe 
that the first 6 functions of Table 4 are the same as those of 
Table 3. The 7™ and 8" functions in [10] were modified 
from those of [7] to properly span the reduced-gradient 
Nedelec space of order 1.5, but while the functions of order 2 
and higher were modified to take on rotational and gradient 
forms, they were not constructed to span the reduced- 
gradient spaces. The edge-based gradient functions (210-412 
and 718-220) are slightly different from those of [7], while 
the face-based functions are substantially different. 


A set of orthogonalized basis functions was proposed in 
2001 by Sun, Lee, and Cendes, and these share many 
functions with the Webb functions of [10]. Table 6 presents 
the functions that differ from those of Table 4 up to degree 3 
(some of the face-based rotational functions). 


Basis function families subsequently proposed by Lee, 
Lee, and Lee [14] and Preissig and Peterson [15] also share 
many members with the Webb bases, and differences 
between the sets of [14] and [15] and those of Table 4 are 
described in Tables 7 and 8, respectively. The functions of 
[14] and [15] differ from the previous sets in that they have 
subsets that properly span both polynomial-complete spaces 


as well as the reduced gradient spaces of Nedelec [1]. An 
analogous set of divergence-conforming vector basis 
functions was proposed for tetrahedral cells by Botha [19]. 
Table 9 reports the Botha functions, after they were projected 
onto triangles and rotated by 90 degrees to obtain curl- 
conforming bases. It is clear that they have much in 
common with the previous families. The Botha functions 
also span both polynomial complete and reduced spaces. 


The functions of [14] and [15] were the first families to 
span polynomial complete and reduced gradient spaces. 
However, condition number studies carried out in [32] show 
that these basis functions exhibit rather poor linear 
independence compared to the basis functions proposed by 
Ingelstróm in 2006 [18], which also span both types of 
spaces. (Although these functions were proposed for 
tetrahedrons, they are easily projected onto triangular 
domains.) The linear independence of the Ingelstróm basis is 
enhanced by ensuring that the projection of the higher-order 
basis functions onto samples of lower-order functions vanish. 
The Ingelstróm bases are presented in Table 10. 


In the 2007 text [20], Demkowicz presented two families 
of basis functions for triangular domains. The family that he 
labels “Nedelec type 1” (reported in Table 11) spans both 
polynomial complete spaces and the reduced gradient spaces. 


Several additional families of hierarchical bases have been 
proposed for triangles that span polynomial complete spaces. 
Cai et al [12] proposed a set of divergence conforming 
functions; those functions have been converted into a curl- 
conforming set and are presented in Table 12. Another set of 
curl-conforming vector basis functions were proposed by 
Ainsworth and Coyle [11], and these are presented in Table 
13. Schóberl and Zaglmayr proposed both curl-conforming 
and  divergence-conforming functions for triangular, 
quadrilateral, tetrahedral, hexahedral, and prism shapes [16- 
17]. In that sense they were the first to cover the realm of 
typical cell shapes in use today. Their curl-conforming 
triangular bases are summarized in Table 14. Table 15 
presents the Demkowicz family [20] labeled “Nedelec type 
2," which properly spans polynomial complete spaces. More 
recently, Xin and Cai [22] proposed a new set of bases 
(Table 16) for polynomial complete spaces. Perhaps as a 
consequence of the fact that these families of functions are 
restricted to polynomial complete spaces, the bases in Tables 
12-16 do not appear to have attracted much interest from the 
computational electromagnetics engineering community. 


Functions that span reduced gradient Nedelec spaces 


The first family of hierarchical basis functions designed 
specifically to span the mixed-order Nedelec spaces was that 
of Andersen and Volakis [8-9]. Reference [8] presented 
functions of up to order 2.5 for tetrahedral cell shapes, while 
reference [9] presented similar functions for triangles. While 
the specific р = 2.5 functions in equations 53, 54, and 55 of 
[9] are linearly dependent and thus invalid, a valid basis for 
triangles can be obtained from the edge-based and face-based 
functions from [8]. Those functions are reported in Table 17. 
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Table 3 
Curl-conforming hierarchical bases on triangles: Wang & Webb 1997 [7] 





Edge Rotational: 


(1) Q, — vVw — wVv 
Q) Q, = wVu — uVw 
(3) Q, 2 uVv - vVu 


complete in reduced gradient space: p — 0.5 








Edge Gradient: (4) V(vw) 
(5) Уби) 
(6) V(uv) 
complete-to-polynomial degree p = 1 
Face (7) uvVw 
(8) vwVu 
(9) wuVv 


Edge Gradient: 


(10) V(uv(u-1/2)) 
о) V(vw(v-1/2)) 
Q2. — V(wu(w-1/2)) 


complete-to-polynomial degree p = 2 








Face 


Edge Gradient: 


(3 uv’ Vw 


(14) wvVw 
(15) vw Vu 
06) wwVu 


(17) WV 


a8 — V(uv(u-1/2)v-1/2)) 
a9 — V(w(v-1/2)(w-1/2)) 
Qo | V(wu(w-1/2)(u-1/2)) 


complete-to-polynomial degree p = 3 
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Table 4 
Curl-conforming hierarchical bases on triangles: Webb 1999 [10] 








Linear bases: 


(1) — (6) are the same functions as in Table 3 











Face Rotational: (7) УО, = uQ, (Q, is defined in Table 3) 
(80 | wQ,—-vQ, 
complete in reduced gradient space: p — 1.5 
Face Gradient: 9) Ушу} 
Edge Gradient: (10) — Vlvw(v-w)] 
aD — ViwuQw-u)] 
42 V{uvu-v)} 
complete-to-polynomial degree p = 2 
Face Rotational: (13) uv(u — v)Vw 
(14) vw(v — w)Vu 
(15) wu(w — u)Vv 
Face Gradient: (16) V{uvw(u — v)) 
(17) V(vwu(v - w)) 
Edge Gradient: (18) V{uv(2u —v)(u — 2v)} 
(19) V(vw(2v - w)(v - 2w)) 
(20) V(wu(2w — u)(w - 2u)) 
complete-to-polynomial degree p = 3 
Table 5 
Webb's orthogonalized functions [10] that differ from those of Table 4: 
i —vw(1398u — 699v + 86)Vu + wu(1398v — 699u + 86)Vv 
S +2020uv(u — v)Vw 
id 2020vw(v — w)Vu — wu(1398v – 699у + 86)Vv 
e +uv(1398w — 699v + 86)Vw 
T vw(1398u — 699w -- 86)Vu 4- 2020wu(w — u)Vv 


—uv(1398w — 699u + 86)Vw 
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Table 6 


Sun, Lee, and Cendes 2001 functions [13] that differ from those of Table 4: 




















(8) wO, (Q, is defined in Table 3) 
(13) vw(v — w)Vu * wu(w — u)Vv + иу(и – у)Уу 
x (393u 4- 80v — 212w)vwVu — (393v 4 80u — 212w)uwV v 
es +(—292u + 292v)uvVw 
i (-131u - 168v — 124w)vwVu 4 (-131v - 168u — 124w)uwV v 
5 
eo) +(-44u — 44v + 262w)uvVw 
Table 7 
Lee, Lee, and Lee’s 2003 functions [14] that differ from those of Table 4: 
Face Rotational: (13) uv(vVw = МУУ) = uv, 
(14) vw(wVu — uVw) - „МО, 
(15) wu(uVv — vVu) = wuQ, 
complete in reduced gradient space: p = 2.5 
Face Gradient: (16) — (17) are the same functions as in Table 4. 
Edge Gradient: (18) — (20) are the same functions as in Table 4. 
complete-to-polynomial degree p = 3 
Face Rotational: (21) uv(u — v.y(vVw — wVv) 
(22) vw(v — w)wVu — uVw) 
(23) wu(w — u)(uVv — vVu) 
(24) uvw(vVw — wVv) 


complete in reduced gradient space: p — 3.5 
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Table 8 
Preissig and Peterson 2004 [15] functions: 





Linear bases: (1) - (6) are the same functions appearing in Table 3 


Face Rotational: (7) u(vVw —wVv)- „О, 
(8) v(wVu – иу) = УО, 


complete in reduced gradient space: p = 1.5 





Face Gradient: (9) is the same function as in Table 4. 
Edge Gradient: (10) — (12) are the same functions in Table 4. 
complete-to-polynomial degree p = 2 





Face Rotational: (13) — (15) are the same functions in Table 7. 


complete in reduced gradient space: p = 2.5 





Face Gradient: (16) — (17) are the same functions as in Table 4. 
Edge Gradient: (18) — (20) are the same functions as in Table 4. 
complete-to-polynomial degree p = 3 








Face Rotational: (21) — (24) are the same functions that appear in Table 7 


complete in reduced gradient space: p = 3.5 








Table 9 
Botha 2007 functions [19], converted to curl-conforming bases on triangular cells: 





Linear bases: (1) — (6) are the same functions as in Table 3, but the equivalent of (4) — (6) on tetrahedrons include a scale 
factor of 3 


complete-to-polynomial degree p = 1 





Face Rotational: (7) — (8) are the same functions as in Table 8 


complete in reduced gradient space: p = 1.5 





Face Gradient: (9) is the same function as in Table 4. 
Edge Gradient: (10) — (12) are the same functions in Table 4. 
complete-to-polynomial degree p = 2 








Face Rotational: (13) — (15) are the same functions in Table 7. 


complete in reduced gradient space: p = 2.5 
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Table 10 


Ingelstróm 2006 curl-conforming hierarchical bases on triangles [18]: 





Linear bases: 


(1) - (6) are the same functions appearing in Table 3 

















Face Rotational: (7) — (8) are the same functions appearing in Table 4. 
complete in reduced gradient space: p = 1.5 
Face Gradient: (9) is the same function as in Table 4. 
Edge Gradient: (10) — (12) are the same functions in Table 4. 
complete-to-polynomial degree p = 2 
Face Rotational: (13) 3vw(v — w)Vu — uw(2v — w)Vv — uv(v — 2w)Vw 
(14) 3wu(w — u)Vv — vu(2w — u)Vw — vw(w — 2u)Vu 
(15) 3uv(u — v)Vw — wv(2u — v)Vu — wu(u — 2v)Vv 
complete in reduced gradient space: p — 2.5 
Face Gradient: (16) — (17) are the same functions as in Table 4. 
. 2 2 
Edge Gradient: a8)  V{uv(u? - 3uv + v?)} 
(19) V] uw – Зу  w?)] 
(20) V{wucw? — 3wu +u’)} 
complete-to-polynomial degree p = 3 
Face Rotational: (21) Suv(u’ — 3uv+v?)Vw—-V {uvw(u? – Зиу+ у? ) 
2 2 2 2 
(22) 5vw(v^ — 3yw + и )Ми – V {uvw(v – Зуу + )} 
2 2 2 2 
(23) 5wu(w^ — 3wu u^ )Vv — V [uvwiw — 3wu +u )} 
(6u — v — wYXv — w)vywVu * (6v — w — u)(w — u)wuVv 
Q4) 


+(6w —u—v)(u—v)uvVw 


complete in reduced gradient space: p = 3.5 
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Table 11 
Demkowicz 2007 hierarchical bases [20], type 1 Nedelec 





Linear: 


(1) – (6) аге the same functions appearing in Table 3. 
complete-to-polynomial degree p = 1 





Face Rotational: 


(7) w 
(8) uQ -vQ 


complete in reduced gradient space: p = 1.5 





Face Gradient: 


Edge Gradient: 


(9) V fuvw} 


(10) VIvw(w-v)] (11) Vfwu(u-w)} 


a2) Vfuv(v = и)! 


complete-to-polynomial degree p = 2 





Face Rotational: 


3 wQ, (14) wu Q, 


5 — wu Q,-wv Q, 


complete in reduced gradient space: p — 2.5 





Face Gradient: 


Edge Gradient: 


(16) Уу? (17) Vfu’vw} 
(18) ve (f= vy ас! 
(19) riu T wy EE 
(20) vi Е uy e 


complete-to-polynomial degree p = 3 








Face Rotational: 


2 


Ю! 


(21) wv 3 
22) my 0. 
Q3 ww Q, 
(24) wu Q, — wy Q, 


complete in reduced gradient space: p — 3.5 
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Table 12 


Cai et al. 2001 hierarchical bases [12], converted to curl-conforming: 














Edge (1) vVw 
(2) —wVv 
(3) wVu 
(4) —uVw 
(5) uVv 
(6) —yVu 
complete-to-polynomial degree p = 1 
Face (7) vwVu 
(8) uwV v 
(9) uvVw 
Edge ao — wViv-u] 
an — wViw-v| 
(12) wuV {u – w} 
complete-to-polynomial degree p = 2 
Face (13) vw(w-— v)Vu 
14) . uw(w—u)Vv 
(15 . uv(v—u)Vw 
(16) uvwV u 
(17) uvwVv 
Edge a8)  uv(v-u)V{v-u} 
(19) vw(w-v)V{w-v} 
(20) wu(u — wv tu — w} 


complete-to-polynomial degree p = 3 
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Table 13 
Ainsworth and Coyle 2001 hierarchical bases on triangles: [11] 





Linear bases: (1) – (6) аге the same functions as in Table 3 but with a scale factor of 2 
complete-to-polynomial degree p = 1 











Face (7) vwVu 

(8) wuVv 

(9) uvVw 
Edge (10) ат! yVw киеу WN V 

(11) ЕЕ ES Sd oue EE 

(12) yy Р Эу, 

complete-to-polynomial degree р = 2 

Face (13) uvwVu 


(14) uvWV v 

(15) | vw(w— v)Vu 
(16) wu(u — w)Vv 
(07) . uv(v—u)Vw 


15(w- vy -4w-v)-5. 15(w- vy - 4(w- v)-5 


Edge (18) z Vw+ UNE wVv 
ad 15(и – и) – 4(и – у) – 5 Wu 15(u— w) +4(u—w)-5 uVw 
6 6 
60) I5(v — u) a u)—5 Vv + I5(v — u) 120 u)—5 Nu 


complete-to-polynomial degree p = 3 
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Table 14 
The Schoberl/Zaglmayr hierarchical basis functions [16-17] 
(scale factors of 2 in many of the original functions are omitted below): 





Linear bases: (1) - (6) are the same functions appearing in Table 3 
complete-to-polynomial degree 1 





Face Rotational (7) w(uVv — vVu) (8) wV (uv) — uvVw 
Face Gradient: (9) is the same function as in Table 4. 
Edge | Gradient: (10) — (12) are the same functions in Table 4. 


complete-to-polynomial degree p = 2 








Face Rotational (13) w(2w — D(uVv — vVu) 
(14) w(2w — DV(uv) - uvV(w(2w — 1)} 
(15) uv(u — v)Vw — wV (uv(u — v)) 


Face Gradient: (16) V{uvw(u — Y) L (17) V{uvw(2w — 1)} 


Edge Gradient: as) V{uv(u? - 3uv + v’)} (19) V{vw(v? — 3vw+ w?)} 


(20) V{wucw? — 3wu + и?) 


complete-to-polynomial degree p = 3 








Table 15 
Demkowicz 2007 hierarchical bases [20], type 2 Nedelec 





Linear bases: (1) – (6) аге the same functions as in Table 3 
complete-to-polynomial degree p = 1 





Face (7) vwVu (8) wuVv 
(9) uvVw 
Edge (10) — (12) are the same functions as in Table 11 


complete-to-polynomial degree p = 2 








Face (13) uvwV u (14) uvwVv 
(45)  wv’Vu (6) — wu Vv 
(47 . uv'Vw 

Edge (18) — (20) are the same functions as in Table 11 


complete-to-polynomial degree p = 3 
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Table 16 
Xin and Cai 2011 hierarchical bases on triangles [22]: 





Linear bases: 


Edge | Gradient: 


(1) - (3) are the same functions appearing in Table 3. 


3 of (wv) 
"i z 2 | 


(5) E ven -w | 





10 2 
(6) У [e] 
2 2 


complete-to-polynomial degree p = 1 





Face 


Edge Gradient: 


(7) eJ5 vwVu 
(8) 6\/5 wuVv 
(9) 3410 uvVw 


(10) £ V(Qw- v» -(w-v)) 
(11) 2 V((u- wy - (u- w)) 
(12) E V((v-uy -(v-u)) 


complete-to-polynomial degree p = 2 








Face Rotational: 


Face 


(13) 4V105 vw(w-v)Vu 
(14) 4V105 wu(u—w)Vv 
(15) 24210 uv(v—- u)Vw 


(16) 124/35 uvwV u 
(17) 12455 uvwV v 
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Table 16, continued 
Xin and Cai [22] hierarchical bases on triangles: 


























N70 5 3 
Edge Gradient: (1800 | -(и-у) - Z(w- vy. 
20 4 2 
470 .( 5 "E я 
(19) —— M —(u- w) ——(u—-w 
20 4 ( ) 2 ( ) 
7 5 4 3 2 
(20) — V/ -—(v-u)' -=(v-u 
ae Vege ) 5 ) 
complete-to-polynomial degree p = 3 
Table 17 
The Andersen and Volakis hierarchical basis functions on triangles [8-9]: 
Edge Rotational: (1) - (3) are the same functions appearing in Table 3. 
complete in reduced gradient space: p = 0.5 
Edge (4) (v—w) Q, 
(5) (и-и) ©, 
(6) (u-v) Q, 
Face Rotational: (7) u О, 
(8) vO, 
complete in reduced gradient space: p = 1.5 
Edge (9) (v-wy Q, 
00 — (w-uy Q, 
a) (u-v? Q, 
Face (12) V(uvw) 


03) w Q, 
(44 vQ, 
as) w° Q, 


complete in reduced gradient space: p = 2.5 
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Table 18 


The Graglia et al. hierarchical basis functions on triangles [21]: 

















Edge Rotational: (1) — (3) are the same functions appearing in Table 3. 
complete in reduced gradient space: p = 0.5 
Edge à — 43(v-w)Q, 
6 — V3(w-u) Q, 
6) — 43(u-v) Q, 
Face Rotational: (7) 243 u О, 
(8) 243 v Q, 
complete in reduced gradient space: p — 1.5 
Edge 6)  S|3v-w?-1-uu-2) Q,/2 
(10) — J5i3w-uy -1-vv-2] Q,/2 
an 45[3Xu-vy -1-wQw-2)] Q,/2 
Face Rotational: (12) 6V5(v — w)u Q, 
(3) 6V5(w-u)v Q, 
04) — 243 u(5u - 3) Q, 
05 — 243 v5v-3) Q, 
complete in reduced gradient space: p — 2.5 
Edge ao HIS WD 23v -w) - 3u(u - 2(v - w)] Q, /2 
an V7 {Sw =u) - 3(w-u) - 3v(v - 20w —u)] Q, /2 
a8 — J7[5(u- v? - 3(u- v) -3w(w - 2u у) Q,/2 
Face Rotational: (19) 24/30 (u^ – 8и + 2)и Q, 
Qo 2430 (1? -8v-2» Q, 
QD — 2430(v-wY7u - 3u Q, 
(22) 24/30(w – и)(7у – 3v О, 
(23) 24210 (? — Avw 4 w^)u Q, 
(24) 24210 (w° — 4wu +u? )v О, 


complete in reduced gradient space: p = 3.5 
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After the publication of [8-9], the basis families proposed 
in [14-15] and [18-20] appeared; these properly span both 
polynomial complete spaces and the reduced gradient spaces 
of [1] and were discussed above. The authors’ own 
contribution [21], which spans the reduced gradient spaces 
(but not polynomial complete spaces) was proposed in 2011. 
Those functions are reported in Table 18. As mentioned 
above, one limitation of hierarchical bases is that they often 
exhibit rather poor linear independence as their order is 
increased. Condition number comparisons carried out in 
reference [32] suggest that the bases of [18] and [21] exhibit 
better linear independence than many of the other families 
that span the reduced gradient spaces on triangles. In 
general, the use of orthogonal polynomials has been found to 
enhance the linear independence of the basis functions. 


We note that a recent paper [24] presents a comprehensive 
review of the construction of various scalar and vector basis 
families, including the triangular cell case. A framework for 
defining curl-conforming bases on triangles is described, 
with enough generality to incorporate various polynomial 
families to enhance linear independence. 


V. QUADRILATERAL BASIS FAMILIES 


Vector basis functions for quadrilateral or hexahedral cells 
were proposed by Ainsworth and Coyle [11], Notaros and his 
colleagues [23-24], Jorgensen et al. [25], Demkowicz [20], 
Graglia and Peterson [28], and Zaglmayr [17]. The functions 
in [24] are divergence-conforming, and in the tables that 
follow these are converted into curl-conforming bases by a 
90 degree rotation of the 2D functions. Tables 19-24 present 
the quadrilateral bases of [11], [23], [24], [25], [20], and 
[28], respectively, on a reference domain -1 € u € 1 and -1 € 
v<1. As in the triangular-cell case, basis functions labeled 
“edge” straddle adjacent cells and must be adjusted to 
maintain tangential-vector continuity across the common 
edge while basis functions labeled “face” are local to a cell. 
For brevity these tables report only the 0 -сотропепі 
functions; the v-component functions can be obtained by 
exchanging u and v in the expressions. 


It is noteworthy that these families all involve polynomial 
representations of the same maximum degree, with 
differences only in the specific form of the polynomial 
function or the scale factors used with the basis functions. 
These basis families all span the same spaces, both 
polynomial-complete and reduced gradient, as indicated in 
Table 19. The differences in specific polynomials and scale 
factors affect the linear independence of the basis functions. 
Reference [20] presented two different families for 
quadrilateral cells, constructed from either Peano or 
Legendre polynomials. The Legendre polynomials are 
believed to provide better linear independence and were used 
to construct the bases in Table 23. Numerical experiments 
[34] show that the functions of [23] exhibit relatively poor 
linear independence, which was improved by the functions 
of [24] and [25]. Of the functions tested, the families in [25] 
and [28] exhibit the best linear independence of the 
quadrilateral group. 


It is interesting that, in contrast to the typical development 
of triangular-cell curl-conforming bases, those in Tables 19- 
24 contain no functions that are explicit gradients. The 
development of basis functions for hexahedral cells that can 
be grouped into “rotational” and “gradient” sub-types was 
proposed in 1998 by Midgard and Sande [35], who only 
presented functions up to degree 2. However, the family 
proposed by Zaglmayr [17] for quadrilateral cells, 
summarized in Table 25, also employs gradients to the 
maximum extent. One benefit of working with explicit 
gradients for the edge-based functions is that, when 
implementing adaptive p-refinement with the vector 
Helmholtz equation, the user has the flexibility of 
introducing individual transition functions into the 
representation without exciting spurious modes. Without the 
use of explicit gradients, the representation used with a cell 
must be complete to the highest order of any of the 
surrounding edges to avoid spurious modes [33-34]. In 
actuality, it is straightforward to convert the edge-based 
functions in the previous tables into gradient form by 
“completing the gradient’ — a procedure similar to 
“completing the square” in algebra — and as an illustration 
Table 26 presents those functions for the family in Table 20. 


VI. APPLICATION TO SURFACE INTEGRAL EQUATIONS 


The functions under consideration should be appropriate 
for use with typical surface integral equations of 
electromagnetics, the EFIE and MFIE. However, research 
suggests that the EFIE only provides a self-consistent surface 
charge density and superconvergent far fields when it is 
discretized with divergence-conforming bases from the 
reduced-curl Nedelec spaces [2-3, 5]. Consequently, 
considering the matrix conditioning comparisons reported in 
[32], the best choices for a triangular-cell basis would appear 
to be the Ingelstróm family [18] or the Graglia et al. family 
[21], suitably converted into a divergence-conforming form. 
For a quadrilateral-cell representation, matrix conditioning 
studies in [33] suggest the use of the Jorgensen et al. bases 
[25] or the Graglia and Peterson [28] bases. Note that the 
bases in [21] and [28] are designed to accommodate both cell 
shapes in a single mesh. 


Recent research also suggests that the optimal basis type 
for the MFIE (to produce superconvergent far fields) should 
span polynomial-complete spaces rather than the reduced- 
gradient or reduced-curl spaces [3-5]. It appears that either 
divergence-conforming or curl-conforming bases may be 
used with the MFIE. Thus, for triangular cell shapes the 
functions of Ingelstróm [18] are suggested (note that the 
Graglia et al. family in [21] does not span polynomial 
complete spaces, and thus may not be optimal for use with 
the MFIE). For quadrilateral-cell representations, both the 
Jorgensen et al. bases [25] and the Graglia and Peterson [28] 
bases do span polynomial complete spaces and should be 
suitable. If it is desired to mix both triangular and 
quadrilateral cell shapes with the MFIE, one of the other 
families may be considered. 
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Table 19 
Ainsworth & Coyle 2001 bases for a quadrilateral domain: [11] 


(Only 4 components reported below; swap u and v to get » -components.) 


























Edge (1) ü (1—v)/2 (2) п (1+у)/2 
4 functions, complete in reduced gradient space: p = 0.5 
Edge (5) ü u(1—v)/2 (6) uu(l+v)/2 
8 functions, complete-to-polynomial degree p = 1 
Face (9) á (v! -1)/2 40  &u(v -D/2 
12 functions, complete in reduced gradient space: p — 1.5 
Edge 43 — & (3 -TD(10-v)/4 14) 0 (3и? –1)(1+у)/4 
16 functions, complete-to-polynomial degree p = 2 
Face 07) & Bu? -1)(v* -1)/4 a8) — &(v —v)/2 
49 — &u(v —v)/2 (20) 0 (3и -D(v-v)/4 
24 functions, complete in reduced gradient space: p = 2.5 
Edge (25)  ü(5w —3uy1—v)/4 Q6)  ü(5w – 3и)(1+у)/4 
28 functions, complete-to-polynomial degree p = 3 
Face (29) й (5и – Зи)(у–1)/4 (30) й (5и —3uyv! —v)/4 
61) 0 (53-62) /8 (32) й и(5у* – 6у>)/8 


зз) 0 (3и? –1)(5у –6у2)/16 (34) 0 (5и – З3и)(5у* – 6у)/16 


40 functions, complete in reduced gradient space: p = 3.5 
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Table 20 


Ilic & Notaros 2003 bases for a quadrilateral domain: [23] (H components only) 








(1) — (12) are the same functions as in Table 19, without the factors of 2 














Edge 43 — &u(-v) 104)  & u (0v) 
Face 00  &uw(v-1) (8) u@ vy) 

(19) H uW" —v) Q0  üu(v -v) 
Edge Q5)  ü (uw —uy1—-v) Q6 й (и`-и)(1+у) 
Face Q9  ü(w-uyv -l1 G0)  & (i —u)yv! —v) 

31)  ü& (v -1 G2)  ü u(v* —1) 

G3  ü(w-DQ'-D G4) й (и –и)(у* – 1) 

ТаЫе 21 
Djordjevic & Notaros bases for a quadrilateral domain: [24] (% components only) 
(1) – (8) аге the same functions as in Table 19, without the factors of 2 

Face (9) ü 2(v —1) a0 — &2u(v 1) 
Edge 43 0 (2и -D(u-v) (14)  & Qu^ - D(ou4- v) 
Face a) a 2(2u* —1)(v’ - 1) (18) 0 dw" v) 

(19) 0 duy" у) со &4Qu -LDOo =) 
Edge Q5)  ü (4u —3u)1- v) Q6 ü (Au – Зи)(1+ у) 
Face Q9  ü2(Aw —3u)(v! - 1) G30) & 4(4u^ — 3u)(v? — v) 

(31) — &(8v' —10v? 42) G2 йи(8у* –10у° +2) 

(33)  & Qu -1(8v -10 -2) (34) й (4и – 3и)(8у* – 10у +2) 
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Table 22 


Jorgensen et al. 2004 bases for a quadrilateral domain: [25] (ù components only) 


(The first 28 functions are the Ainsworth & Coyle bases with different scale factors.) 








Edge 


Edge 


Face 


Edge 


Face 


Edge 


Face 


~ 13 
1 —(1- 
(1) 1 v) 
(5) U ED 
a 588 
9 — (v -1 
(9) u 255 ) 


15 55 
i -1)ü-» 


= 


(13) 


= 


(17) БО —1)(v? -1) 


a9 0 JZ о) 


„121 
Q5 ù ig (5 -3u)(l-v) 
. [105 
Q9 й is P" —3и)(°° 1) 
GD) й si; 6v -6v «1 


A. 13 
2 — (1+ 
(2) u 435 v) 
(6) U I] 
(10) a Euo –1) 


15 


04 й ig ae» 
(8) 0 (и) 

оу | >(3u?-1)(v*-») 
Q6 ù n oue) 
(30) ù 7 (5u’—3u)(v" =») 


2 1945 4 3 
(32) u 315 0 — бу +1) 


63 ù үз (зи -I)(5v* - 6v? 1) 


„ |2205 


G4 — f J—— (5w -3u)(5v* - 6v? 1) 


2048 
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Table 23 


Demkowicz 2007 bases for a quadrilateral domain: [20] (7% components only) 


(The first 28 functions are the Ainsworth & Coyle bases with slightly different scale factors.) 








(1) — (16) are the same functions as in Table 19. 








n ài -1(v^-1 . Q- 
Face (17) u оо (18) u е0 
5 4 4 
3 ЕК С 
(19) ü jee (20) ü 3u yr) 
4 2 4 
Edge (25) — (28) are the same functions as in Table 19. 
„ 5и? — 3u (v? —1) „ 5и? – Зи (у - v) 
Face (29) EN (30) A ——— ——— 
2 4 2 4 
.(?-D(5v -1 2-2] 2—1 
sy. еа д, Bye Or sD 
4 4 4 4 
, à -1(? -1) (5? -1 
(33) Р и (v ) (5v ) 
2 4 4 
, 9v -3u (v! - D) (5 -1 
(34) Р и и (у ) (Sv ) 
7 4 4 
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Table 24 


The Graglia & Peterson hierarchical basis functions [28]: 


(Only 4 components reported below; swap wu and v to get » -components.) 








Edge 


Edge 


Face 


Edge 


Face 


Edge 


Face 


(0 H (V) 
(5) й V3u(1—v) 


Jio 


10 
б, a-v 
(9) и 5 v“) 


aj 





2 
43) ù E а-ы 





a) ù Ler -Dü-») 
49) ù 2m uv(1— v?) 


Q5  &-(5w -3u)(1-v) 


IPS 


Q9 ù (5^ - 3u)(v? -1) 


su 
© 


A 


в) й J[— (7v! - 8v «1) 


E 
[ul 


Q) ü (19v) 
(6) й V3u(1+v) 


. J30 


10 Sc A my? 
(10) "та. v^) 





a4) ù 5|® Hew 


JO 


2 





a8 ù v(1—»?) 


ч 


оой E —1)у(1—›?) 


(26)  &—-(5w -3u)(1ev) 


FE 


, N490 


(30) u ME -3u)(v? -v) 


G2 ù [unici +1) 


Q3 й [Be i-e) 


„ [105 


(34) ù — (5i -3u)(7v* - 8v? +1) 


52 
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Table 25 
The Schoberl/Zaglmayr 2006 bases for a quadrilateral domain: [17] 


(Roughly half the functions are shown of the first 24; swap u and v, à and Y to get their pairs, except as noted.) 








Edge: | (D à (1—v) Q) ü (1v) 
Edge gradients: 

5 &ul-v)-? 0-w)/2 (6) ü u(1-v)-$ —w)/2 
Face: (9) | & (v -1D 40 H uW" —-D-$ (u^ -Dv 


u(v? —1)+0 U - Dv 


= 


The pair for (10) is the gradient: (12) 


Edge gradients: 
(13) #0 (1—-v)(3u? - D-$ Qi? —u) (14) 0 (1+У)(3и2-1)+7 (и-и) 


Face (7)  ü& (3) -D^ —D/24$ QP? —u)v 

a8) й (У-у) 

(19) йа(у -v)-?^-DGv -D/2 
(34^ е1) (у —-yv)-$9(Qv -DüG£ —u) 


= 


(20) 


The pair for (20) is the gradient: (24) a (3u’—1)(v’—v)+0 (3v? —1)(u° -u) 








Table 26 
Ilic & Notaros edge basis functions converted to gradients. These functions may be used to replace those in Table 20. 








(5) й 2и(1-у)+? (1-и?) (6) й 2и(1+у)-? (1-и?) 
(13) 0 (3и -1)(1-у)-–? (и-и) (14) й (3и –1)(1+у)+7 (и? –и) 
(05) й 4(и`-и)(1-у)-– (и – 2и?) (06) й 4(и`-и)(1+у)+? (и – 2и?) 
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VII. CONCLUSIONS 


This article has attempted to organize and review the 
variety of different types of hierarchical vector bases 
appropriate for use in solving the common equations of 
electromagnetic fields on (possibly curved) triangular or 
quadrilateral domains. Tables are used to illustrate 16 
different basis families for triangular domains and 7 families 
for quadrilateral domains. There is considerable variety in 
the specific form of the triangle basis functions, while the 
quadrilateral basis families are much more similar to each 
other. Some relative advantages are identified that may offer 
benefits for use with the numerical solutions of surface 
integral equations. 
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